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Rocking and rolling: a can that appears to rock might actually roll 
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A beer bottle or soda can on a table, when slightly tipped and released, falls to an upright position 
and then rocks up to a somewhat opposite tilt. Superficially this rocking motion involves a collision 
when the flat circular base of the container slaps the table before rocking up to the opposite tilt. 
A keen eye notices that the after-slap rising tilt is not generally just diametrically opposite the 
initial tilt but is veered to one side or the other. Cushman and Duistermaat (2006) recently noticed 
such veering when a flat disk with rolling boundary conditions is dropped nearly flat. Here, we 
generalize these rolling disk results to arbitrary axi-symmetric bodies and to frictionless sliding. 
More specifically, we study motions that almost but do not quite involve a face-down collision of the 
round container's bottom with the table-top. These motions involve a sudden rapid motion of the 
contact point around the circular base. Surprisingly, like for the rolling disk, the net angle of motion 
of this contact point is nearly independent of initial conditions. This angle of turn depends simply 
on the geometry and mass distribution but not on the moment of inertia about the symmetry axis. 
We derive simple asymptotic formulas for this "angle of turn" of the contact point and check the 
result with numerics and with simple experiments. For tall containers (height much bigger than 
radius) the angle of turn is just over n and the sudden rolling motion superficially appears as a nearly 
symmetric collision leading to leaning on an almost diametrically opposite point on the bottom rim. 

PACS numbers: 45.20.dc,45.40.-f 
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I. INTRODUCTION 



After a meal or a drink, fidgety people sometimes play 
with the available props, such as empty glasses, bot- 
tles, and soda cans. These containers are all nearly axi- 
symmetric objects with round bottoms. A natural form 
of play is to roll these containers on their circular bot- 
toms. When a slightly tipped container (Fig.[T]) is let go, 
sometimes one sees it fall upright, make a slight banging 
sound, and then tip up again at another angle, then fall 
back upright again, and so on for a few rocking oscilla- 
tions. Because these oscillations usually damp quickly its 
easy to miss the details of the motions. To aid the eye, 
instead of just letting go of the slightly tilted container, 
you can flick its top forward with the fingers. This pro- 
vides an initial righting angular velocity along the axis 
about which the container was initially tipped (too large 
an angular velocity will cause the container to lift off the 
table as it pivots over). The tipping container, as before, 
falls to a vertical configuration, at which time its bot- 
tom circular face bangs on the table. After this bang, 



the container tips up onto the other side and maybe falls 
over. When this experiment is performed with a con- 
tainer that is not too tall and too thin, you can see that 
the container does not fall exactly onto the diametrically 
opposite side of the bottom rim. That is, point A on the 
bottom of the container that initially contacted the table 
and the new contact point B, that the container rocks 
up onto, are not exactly 180 degrees apart. This exper- 
iment (video available online Q) is shown schematically 
in Fig. [2 Fig. [3] shows a histogram of one particular con- 
tainer's orientations after we repeatedly flicked it. The 
distribution is strongly bimodal with no symmetric falls 
over many repeated trials. 
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Note the apparent symmetry breaking with Ayb ^= n. 
Is this deviation from symmetric rocking due to imper- 
fect hand release? Here we show that the breaking of 
apparent symmetry is consistent with the simplest deter- 
ministic theories, namely smooth rigid body dynamics. 
We derive formulas for the "angle of turn" , for both the 
perfect-rolling and for the frictionless-sliding cases given 
infinitesimal symmetry breaking in the initial conditions. 
The results here generalize some of the nearly-falling-fiat 
results of Cushman and Duistermaat Q ; they considered 
the special case of the pure-rolling of flat disks. 
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FIG. 1: Releasing a tipped container. After this release 
the container falls and rocks as shown in Figs. [2] and [3] 
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FIG. 2: Toppling a container. Schematic time-sequence of 
container (illustrated as a cylinder) during one tipping exper- 
iment: side view (top row) and top view (2nd row), a) Tilt a 
container and flick its top with your finger, so that b) the con- 
tainer becomes vertical again, and then c), d) tips part way 
over some more and then eventually d) has sufficient energy 
to fall on its side. The container hardly ever falls symmet- 
rically to the diametrically opposite side with Alp = n. See 
Fig! 



II. CANDIDATE THEORIES FOR SYMMETRY 
BREAKING 



If the table was perfectly planar and horizontal, the 
container's bottom was perfectly circular, and the con- 
tainer was perfectly axisymmetric, then an initial condi- 
tion with a purely righting angular velocity would re- 
sult in a collision in which, just as the container be- 
comes vertical, all points on the container's bottom slap 
the table-top simultaneously. The consequence of such a 
rigid-container collision is not computable since algebraic 
rigid-body collision laws are not well-defined for simulta- 
neous multi-point collisions; for example,the order of the 
impulse locations is then ill-defined [l|, 0, H, [H . Basing 
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FIG. 3: Typical toppled container orientations. A his- 
togram of the directions along which a container fell for 42 
trials. The container always falls to one side or another, 
but essentially never to the diametrically opposite side cor- 
responding to 180 degrees. The distribution about Atp = 180 
degrees is strongly bimodal. Whether the fall is to the left 
or to the right depends sensitively on initial conditions. The 
leftward falls have a mean of about Atp — 180 + 33 degrees 
and the rightward falls have a mean of about Aip = 180 — 37.9 
degrees. For this cylinder theory predicts Aip = 180 ± 40 de- 
grees for frictionless sliding and Aip — 180 ± 22 degrees for 
rolling with no slip. 



the rocking outcome on such a perfect flat collision would 
be basing the outcome on details of the deformation, that 
is, leaving the world of rigid-object mechanics. 

Geometric perturbations to the circular rim can sim- 
ilarly lead to two-point collisions, where both points of 
contact are now on the rim. The second point of con- 
tact could be diametrically opposite to the first point, or 
anywhere else on the rim. We could then compute the 
consequences of, say, a plastic collision at the new point 
of contact. However, the consequences of the collision 
would depend on the location of the geometric imper- 
fection. For a given collision point, such a theory could 
predict the energy dissipation at the collision. But such a 
theory cannot be useful in predicting a systematic break- 
ing of symmetry as seen in Fig. [3] For any given imper- 
fection, the motion is deterministic and depends on the 
location of the impection, and there is no reason to ex- 
pect that the location of the imperfection would have a 
distribution similar to that in Fig. [3] 

Assuming a perfectly flat rigid bottom, the circle- 
slapping-ground simultaneous collision is essentially im- 
possible. After all, the container is being launched by 
imperfect human hands that cannot provide any exact 
initial conditions. Accounting for various symmetries in 
the problem, the set of all motions of a container rolling 
without slip is three dimensional; the space of solutions 
could be parameterized by, say, the minimum tip angle 
4> m in, the yaw rate ip a t that position, and the rolling 
rate 9 at that position. The set of solutions that leads 
to a face-down collision can be characterized with only 
two parameters, a set with co-dimension one 0, 0, 0). 
So, small generic perturbations of a "collisional" initial 
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condition result, with probability 1, in a non-collisional 
motion described by the smooth dynamics of the con- 
tainer rolling or sliding on its circular bottom rim. 

The rest of this paper is about the near-collisional mo- 
tions. We will see that these near-collisional motions 
involve rapid rolling or sliding of the container on its 
bottom rim, so that the contact point appears to have 
switched by a finite angle that is greater than 180 de- 
grees. For this analysis, we assume a geometrically per- 
fect container (axisymmetric) and table (flat). 



III. DYNAMICS OF ROCKING BY ROLLING 

Consider a container with mass m, bottom radius R 
[13], and the center of mass at a height H from the bot- 
tom. The moment of inertia is C about it's symmetry 
axis and is A about any axis passing through the center 
of mass and perpendicular to the symmetry axis. For the 
disk of 0, H = and A = C/2. In our case, H > and 
A > C/2. 

The center of mass of the container is at (xcUGtZg) 
in an inertial frame e x -e y -e z at rest with respect to the 
table. The reference orientation of the container is ver- 
tical as shown in Fig. Any other orientation of the 
container can be obtained from the reference orientation 
by a sequence of three rotations, defining corresponding 
Euler angles as illustrated in Fig. 2J The first rotation 
('yaw' or 'steer') is about e z by an angle ip. This rotation 
also transforms the inertial coordinate axes to e x \-e v \- 
b z \. The second rotation ('pitch' or 'tilt') by an angle 
—<p about the e y i axis results in the e X 2-e y 2-e Z 2 frame 
and determines the orientation of the container up to a 
rotation 9 about the body-fixed symmetry axis (e Z 2 axis) . 
The angular velocity of the container in the rotating e X 2- 
Vyi-^zi frame is therefore entirely along the symmetry 
axis e Z 2 and this relative angular velocity magnitude is 
denoted by 9. 

We will mostly consider two simple extremes for the 
frictional interaction between the table and the con- 
tainer's bottom, namely, sliding without friction and 
rolling without slip. 

For pure rolling, the Euler angles and their first and 
second time derivatives determine the center-of-mass po- 
sition (relative to the contact point), the center-of-mass 
velocity and acceleration, as well as the container's angu- 
lar velocity and angular acceleration. The three second 
order ODEs that determine the evolution of the orienta- 
tion ip, 4>, and 9 follow from angular momentum balance 
about the contact point: 

Qni> + Qi2$ + Qi30 = Si, 1 = 1,2,3. (1) 




e x1 



FIG. 4: Definition of the Euler angles and coordinate axes 
used to define the orientation of the container. 

where 

Qu = A sin 4> — mH R cos <fi + mH 2 sin <p, 

Q12 = 0, Qi 3 = -mHR, 

Qai = Q23 = 0, Q22 = -mR 2 - mH 2 - A, 

Q31 = C cos <fr + mR 2 cos (p — mRH sin (/>, 

Q32 = 0, Q 33 = C + mR 2 , 

51 = (C -2A-2m,H 2 )^coscp + C^9 (2) 
— 2mH Rip<p sin <f> 

5 2 = (C - A + mR 2 -mH 2 )ip 2 sin 20/2 
+ (C + mR 2 )94) sin + mHRip9 cos 
+ mH Rijj 2 cos 2(f> 

+ mg(R cos <f> — H sin 0), 

53 = Cip<P sin (f> + 2m Rip ip(R sin <p + H cos^). 

These equations are derived in Appendix [Al Because the 
center-of-mass velocity is determined by the Euler angles 
and their rates, once these are known the center of mass 
position can be found by integration. 

For the frictionless sliding of the cylinder on its circular 
bottom, the equations of motion are also given by Eq. Q] 
but with 

Qu = Asin0, Q12 = Q13 = Q21 = 0, 

Q22 = -A + mHR sin 2</> - mH 2 sin 2 4> - mR 2 cos 2 0, 

Q23 = 0, Q31 = C cos0, Q32 = 0, Q33 = C, 

51 = -2A^cas<j) + C<j)6 + C(l4cos<l), (3) 

5 2 = (1/2) (C- A)^ 2 sin2</) 

+ C9ip sin <f> + mgR cos <j) — mgH sin <j> 

- mHRtf cos 2<j> + (l/2)m(H 2 - R 2 )<j) 2 sin 20, 

53 = Cip<p sin <(>• 

When the table is frictionless, the horizontal velocity 
of the center of mass is a constant, so that the horizontal 
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FIG. 5: Results of a typical simulation of a container's face 
almost but not quite falling flat on the table-top. Physical 
parameters correspond to the can shown in Fig. \T\ and used 
for the experiments in Fig. [3] m = 0.070; A = 5.1x 10~ 3 ; C = 
8 x 10~ 3 ; g = 9.8; R = 5.1 x 10~ 2 ; H = 6.9 x 10" 2 in consistent 
SI units. Initial conditions ip(0) = 0,-0(0) = 0.005,0(0) = 
tt/6, 0(0) = -0.5, 0(0) = 0, 8(0) = 0.005. All the angles are in 
radians. The step-change in ip is the angle of turn A-0, and 
is a little more than ir as noted in the caption for Fig. [3] 



position is independent of the orientation. The vertical 
position of the center of mass is simply given by zq — 
H cos cp + Rsincp. Without loss of generality the center- 
of-mass can be taken as on a fixed vertical line. 

Our initial discovery of the phenomenon described 
here was found in numerical simulations of the equations 
above. We used an adaptive time-step, stiff integrator 
because the solutions of interest have vastly changing 
time-scales. The singularity of the Euler-angle descrip- 
tion in the near-vertical configuration also contributes to 
the stiffness of the equations. 

For simplicity of presentation, we first describe in de- 
tail the results of integration when the container does 
not slip with respect to the table. First, we integrate 
the equations with initial conditions that lead exactly to 
a face-down collision ip(0) = 9(0) = 0, < 0, and say, 
< < ir/2. As would be expected, the container lands 



in a manner that all of its bottom face simultaneously 
reaches the horizontal table-top. Up to this time the 
motion of the container is identical to that of an inverted 
pendulum hinged at the contact point A. The integration 
becomes physically meaningless at that collision point. 

In the vicinity of these collisional motions, the solu- 
tions of the equations of motion here do not have smooth 
dependence of solutions on initial conditions. To obtain a 
near-collisional motion, we set ip(0) = 0(e), 9(0) = O(e), 
4> < 0, and < 0(0) < 7r/2, where e is a small quantity. 
The results of the integration are shown in Fig. The 
plot of 4>(t) suggests a motion in which the container's 
bottom face periodically comes close to touching the ta- 
ble (0 — > 0), but then gets "repelled" by the floor as if 
by an elastic collision, so that the container rocks down 
and up periodically, ad infinitum, without losing any en- 
ergy, as expected from this dissipation-free system. We 
notice that when « 0, changes almost discontinu- 
ously. Also, when is close to zero, both ip and 9 blow 
up to very large values, resulting in almost discontinuous 
changes in the corresponding angles ip and 9. 

Note that we are simply simulating the apparently 
smooth differential equations, and not applying an alge- 
braic transition rule for a collision. As e goes to zero, the 
angle rates ip and 9 grow without bound, but the magni- 
tude of the angular velocity vector is always bounded, as 
it must be since energy remains a constant throughout 
the motion. In particular, while the angle rates 9 and ip 
are large, they are very close to being equal and opposite 
(0« -V>,Fig.[S). 

Let us examine the consequences of a rapid finite 
change in ip when w 0. The position of the contact 
point Y(xp,yp 1 zp) on the ground relative to the center 
of mass (xg, Vg, zg) is given by the following equations: 

xp — xg — R cos cos ip + H sin 0, 

VP = Vg - Rsinip, 

zp — zg — i?sin0cos"0 — if cos = 0. 

When w and ip is not particularly close to any mul- 
tiple of tt/2, the contact point position is given by 



xp w xg — R cos ip 
VP ~ Vg - R sin ip 
zp pa zg — H = 



(4) 



Given that xg, Vd z g do not change much during the 
brief near-collisional phase (because center of mass ve- 
locity is finite), we can see from Eq. 3] that a rapid con- 
tinuous change in ip corresponds to a rapid continuous 
change in the contact point in a circle with the center 
0(xg, VG, 0). Thus the "angle of turn" AOB defined ear- 
lier is simply the change in ip. At the singular limit of 
a near-collisional motion arbitrarily close to a collisional 
motion, this continuous but steep change in ip approaches 
a step change - this is the "limiting angle of turn" and 
we denote this by Aip tU rn- 
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We consider the near-collisional motions of the no-slip 
container that can be characterized as being the pasting- 
together of two qualitatively distinct motions of vastly 
different time-scales: 

1. inverted pendulum- like motion about an essentially 
fixed contact point when the tilt angle <p is large. 

2. rapid rolling of the container which accomplishes 
in infinitesimal time, a finite change in the contact 
point, and a sign-change in the tilt rate <j). 

Note that the near-collisional motion for a container 
rolling without slip requires very high friction forces (for 
non-zero H, see Fig. [S^). However, plotting the ratio of 
the required friction forces with the normal reaction, we 
find that only a finite coefficient of friction is required for 
preventing slip even in the collisional limit (Fig. [Sfr) . 

The other extreme of exactly zero friction is similar. 
Here, the horizontal component of the center of mass 
velocity may be taken to be zero. The near collisional 
motions for a container sliding without friction are again 
characterized as consisting of two qualitatively different 
phases: 

1. a tipping phase when (f> is not too small, involv- 
ing the container moving nearly in a vertical plane, 
the center of mass moving only vertically, and the 
contact point slipping without friction. 

2. a rapid sliding phase in which the sign of cj> is 
reversed almost discontinuously, and the contact 
point moves by a finite angle in infinitesimal time. 

The angle of turn for the frictionless case and the no- 
slip case are different in general (when H ^ 0). In the 
next two sections, we derive the formulas for the angle 
of turn by taking into account the two-phase structure of 
the near-collisional motion. 



with vastly different time-scales, each regime simple to 
analyze by itself. The overall motion can be obtained 
approximately by suitably pasting together the small-0 
and the large-0 solutions. But the angle of turn is en- 
tirely determined by the small-</> regime, as will be seen 
below. 

First, consider Eq.[T]with i = 1 in the limit of small </>, 
so that we can use sin <j> = <f> and cos 0=1, and generally 
neglect terms of O ((/)). We obtain 

(A(f> - mHR + mH 2 cj))i) - mHRO = 

(C -2A-2mH 2 )ip<j> + Cd(l> (5) 

Now, considering Eq. [T]with i = 3 in the limit of small 
</>, we obtain 

(C + mR 2 - mRH(f>)i> + (C + mR 2 )9 = 2mRH^. 

(6) 

Eqs.[5]and[S]are linear and homogeneous in tp, 9, and 
their time derivatives. So, positing a linear relation be- 
tween ip and 9, we find that the following two equations 
are equivalent to Eqs. [5] and [6) 



= o 



and 



t- = -2- 



(7) 



(8) 



Note that Eq. [7] agrees with the results of the numerical 
simulation, as in Fig. [5j. Even though this equation was 
derived for small <p, this equation is approximately true 
at large (j> as well if the initial conditions of the motion 
at large <j> have ip(0) = O(e) and 0(0) = 0(e), as is the 
case for the near-collisional motions we consider. 
Integrating Eq. [8j we obtain 



(9) 



IV. ANGLE OF TURN FOR NO-SLIP ROLLING 

The rigid body dynamics of disks, containers, and 
similar objects with special symmetries, have been dis- 
cussed at length by a number of authors, including dis- 
tinguished mechanicians such as Chaplygin, Appell, and 
Korteweg. Their works include complete analytical char- 
acterizations of the solutions to the relevant equations 
of motion, typically involving non-elementary functions 
such as the hyper-geometric. Reasonable reviews of such 
literature can be found, for instance, in [5| and [g|. We 
will not use these somewhat cumbersome general solu- 
tions but will analyze only the special near-collisional 
motion of interest to us. 

The calculation below may be called, variously, a 
boundary layer calculation, a matched asymptotics calcu- 
lation (but we are not interested in an explicit matching) 
or a singular perturbation calculation. Essentially, we 
take advantage of the presence of two dynamical regimes 



Thus, when <fi <C 1, both ip and 9 (= —ip) become very 
large. 

Now consider Eq. [1] with i = 2 and with <p -C 1 . 

- (A + mR 2 + mH 2 )4> = ip 2 (l)(C -A + mR 2 - mH 2 ) 
+ ip6<j)(C + mR 2 ) + mHRip9 + mHR-ip 2 + mgR (10) 

Using Eq.[7]in Eq. 1101 and ignoring higher order terms in 
<j), we obtain 

[A + mR 2 + mH 2 )4> = (A + mH 2 )cf)ip 2 - mgR (11) 

Using Eq. [5] for tp in the above equation and neglecting 
the mgR term (as it is a constant and therefore, much 
smaller than l/4> 3 for small <j>), we get 

2 A + mH 2 b\ 



A + mR 2 + mH 2 (j) 3 

h 

/A3 ' 



(12) 
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where 



A + mH 2 , 9 
lf z 

A + mR 2 + mH 2 2 ' 



(13) 



The general solution for the differential equation in 
Eq. [TJ1 can be written as 



4> 2 = rt + b 3 (t-t c ) 2 /0 2 c 



(14) 



where <p c is the lowest (f> attained by the container before 
rising back up again, and t = t c is when this minimum <p 
is attained. Substituting this equation in Eq.[9l we obtain 
a simple equation for the evolution of ip when near the 
surface: 



A2 



<p> 2 <t% + b 3 {t-t c yi<% 



We can now compute the change in ip during a small time 
interval At centered at t — t r . 



Aipt 



-At 

t c +At 
t c -At 

t c +At 

b 2 



ipdt 



tan 



1 {VhAt\ 



dt 



A + mR 2 



iH 2 



A + mH 2 



tan 



/boAt 



From conservation of energy, we can show that b 2 , and 
therefore y/b~3, must be 0((f> c ) (see Appendix IB"]) . Hence, 
VbH/(f>c — > 00 as 4>c — * 0. Keeping At a small constant 
as we let <fi c — > 0, thus approaching the collisional limit, 
gives us the following expression for Aip. 



Aipt 



A 



iH 2 



iR 2 



iH 2 



tan 1 (oo) 



A + mH 2 



A 



iR 2 + mH 2 



(15) 



This is the limiting angle of turn when the container rocks 
by rolling without slip. Earlier numerical integrations 
agree quite well with this formula near the collisional 
limit, as they should. 



Here, bi = 0(e) <C 1 because the initial conditions sat- 
isfy 0(0) = 0(e) and ^>(0) = 0(e) by assumption, as 
before. Using Eq. [16] in the frictionless equation Eq. [T] 
and Eq. [3] corresponding to i = 1, we have, after some 
simplifications: 



t = - 2 i 

ip <p 



(17) 



where b 2 is a constant of integration, whose order is es- 
timated in Appendix [51 Substituting this into the i = 2 
frictionless equation and simplifying by neglecting all 
higher order <p terms, we eventually obtain 



b 2 



A + mR 2 



A3 • 



(18) 



Using arguments identical to the pure-rolling case, this 
results in the following expression for the angle of turn, 



Atpturn = 7T 



for the frictionless limit. 



iR 2 



A 



(19) 



VI. AN ALTERNATE HEURISTIC SMALL 
ANGLE TREATMENT 

In this section, we derive the same angle of turn formu- 
las fEq. [T51 and [H))) without referring back to the compli- 
cated full dynamical equations. Rather, using heuristic 
reasoning, we directly derive equations of motion that 
apply at the small angle limit (0 < <p -C 1). 

We represent the lean of the cylinder by the vector 
$ with magnitude equal to <p> and direction along the 
e y i axis: $ = (pe v \. These quantities and ip are 
related to each other exactly like r, r and 0, respectively, 
in traditional polar coordinates. 

Neglecting any angular velocity component along e z , 
the angular velocity vector for the cylinder is given by 
$ = <pe y i — cjnp&xx. The rate of change of angular mo- 
mentum He about the center of mass G is given by: 



V. ANGLE OF TURN FOR FRICTIONLESS 
SLIDING 

We now briefly outline the procedure for deriving the 
angle of turn for the frictionless case. The procedure 
closely parallels that described for the no-slip case, except 
for small differences below. The frictionless equations 
Eq. [TJ and Eq. [3] corresponding to i = 3 simplifies to 



H G = A<5> = A— (0e w i) 

= -A{0iP + 2<Pi]j)e xl +A{4>-^ 2 )e yl . (20) 

The angular momentum balance equation is then given 
by 



H G = M G , 



(21) 



C—(ipcos<j) + 6) = 
+ ip cos <p = constant = b\ , say. 



(16) 



where M G is the moment of all the external forces about 
G. M G depends on whether or not there is friction; so 
we treat the two cases in turn. 
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A. Sliding without friction 

The vertical position of the center of mass is zq = 
Rs'm(c/)) + Hcos(cj)). So zq w Rip, The vertical ground 
reaction is thus mzg — mRcf), neglecting gravity mg 
in comparison. The moment Mg of this vertical force 
about the center of mass G is equal to —m(Rcos(f> — 
H sixi(j))R(f) « —MR 2 4> in the direction e yl . Substituting 
this along with Eq.f^Olin the angular momentum balance 
equation Eq. 12 1 1 we have 

- MR 2 <f>e y i = A(<f> ~ (j)ip 2 )e y i - A(4>4> + 2#)e xi . 

This vector equation is identical to Eqs. Q~7] and [TH] and 
therefore, lead to the same angle of turn fEq. IT5|) . 



of 202 degrees for rolling without slip. These angles of 
turn would manifest as a deviation of either 40 degrees 
or 22 degrees from falling over to exactly the diametri- 
cally opposite side. In the toppling experiment of Fig. O 
the container orientation was 33 degrees on average from 
180 degrees in the leftward falls (suggesting an angle of 
turn of 213 degrees) and was about 37.9 degrees on av- 
erage from 180 degrees in the rightward falls (suggesting 
angle of turn of 217.9 degrees). The standard deviations 
were respectively 3.9 and 4.5 degrees respectively. The 
experimental angle of turn seems better predicted by the 
asymptotic formula for frictionless sliding in this case. 
Although, neither the frictionless limit nor the no-slip 
limit is just right, both limits capture the many qualita- 
tive aspects of the motion quite well. 



B. Rolling without slip 

The position of the center of mass G with respect to 
the point P' on the cylinder in contact with ground is 
given by r G w r P , - (R - H(j))e xl + (R(f> + H)e z . The 
velocity of the center of mass G is given by r G — r P' + 
$ x (r G — rp/). Using the no-slip constraint Ypi = 0, we 
obtain after some simplifications, the acceleration of the 
center of mass to first order to be 

r G = H(4> - H 2 )e xl + H{cj4 + 2^p)e yl + R4>e z . (22) 

The ground reaction force, which now includes a hori- 
zontal friction force as well, is simply mig, neglecting 
gravity. The moment of this ground reaction force about 
G is given by 

M G = (rpi — r G ) x mrg 

« -B xl H 2 {<j4 + 244) ~ e y i(R 2 4> + H 2 {4> - <j)^ 2 )) 
-e x iRH((f>i) + 2H>). 

Equating M G above to H G from Eq. [20] gives Eqs. [5J 
and [TT]from the previous version of the derivation, there- 
fore resulting in the same formula for the angle of turn 
(Eq.EJ. 

VII. QUANTITATIVE COMPARISONS 

Firsly and most significantly, note that the limiting 
angle of turn does not depend on the initial conditions 
such as the initial tilt and tip velocity. This means that 
we do not have to control these accurately in an experi- 
ment. This also agrees with the relatively small variance 
in the histogram of Fig. [31 in which we did not control 
the initial conditions. 

The histogram Fig. [3] was obtained using the cylindri- 
cal container shown in Fig. [T] with R = 5.1 cm, H = 6.9 
cm, A/m — 5.13 x 10~ 3 m 2 . Using these numbers in the 
angle of turn formulas gives an angle of turn of about 
220 degrees for frictionless sliding and an angle of turn 



VIII. SPECIAL LIMIT: TALL THIN 
CONTAINERS 

For tall thin cylinders with A ~ mH 2 and H >> R, 
both equations for the angle of turn, Eq. [T5] and Eq. 
tend to 7r radians. That is, very tall cylinders are pre- 
dicted to have a smaller symmetry-breaking. This pre- 
diction agrees with the common experience that when we 
tip a tall-enough cylinder (something that looks more like 
a tall thin beer bottle) in a manner that its bottom sur- 
face nearly falls flat, the cylinder essentially rocks up on a 
contact point almost diametrically opposite to the initial 
contact point. Thus, this apparently almost-symmetric 
rocking is accomplished by a rapid asymmetric rolling 
or sliding of the container over roughly one half of its 
bottom rim! This result is basically independent of the 
container-table frictional properties. 

IX. SPECIAL LIMIT: DISKS 

We may define a disk as a container with zero height: 
H = 0, a good approximation is the Euler's disk 
If the radius of gyration of a disk is R g — kR, then 
C = mR 2 k 2 and A — mR 2 k 2 /2. Substituting these into 
the angle of turn formulas for the no-slip (Eq. [T5)) or the 
frictionless (Eq. [T^)) cases, we obtain the same angle of 
turn: 

AVw = irV2 + k 2 /k. (23) 

Indeed, numerical exploration with Eqs. [l][3l suitably 
modified for frictional slip and specialized to disks, shows 
no dependence of the angle of turn Aiptum on the form 
of the friction law or the magnitude of the friction. This 
lack of dependence on friction could be anticipated from 
our small angle calculations for pure rolling cylinders, 
specialized to disks. In particular, we find that the accel- 
eration of the center of mass for a pure rolling cylinder 
(Eq. |22|) with H = is vertical in the small-angle limit, 
enabling the no-slip condition to be satisfied even with- 
out friction. Thus, the rolling solution is obtained with 
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or without friction. For the special case of pure-rolling 
disks, Eq. [221 was found in Q. 

A homogeneous disk has k = l/y2. The corre- 
sponding angle of turn is equal to ny/E « 2.23-7T, which is 
about 41 degrees more than a full rotation of the contact 
point. This prediction is easily confirmed in casual ex- 
perimentation with metal caps of large-mouthed bottles 
or jars on sturdy tables - for such caps, we observe that 
the new contact point is invariably quite close to the old 
contact point. 

A ring such as the rim of a bicycle wheel has k w 1. 
The corresponding angle of turn is equal to ny/3 « 1.737T, 
which is about 48 degrees less than a full rotation of the 
contact point. Thus the apparent near-collisional behav- 
ior of a homogeneous disk and a ring will be superficially 
similar, even though the actual angles of turn differ by 
about 90 degrees. 

The angle of turn can be controlled by adjusting k. For 
instance, between a ring and a disk is an object that ap- 
pears to bounce straight back up the way it falls. And it 
is possible to increase the theoretical angle of turn with- 
out bound by choosing k — > 0, a disk in which almost all 
the mass is concentrated at the center. 



X. SUMMARY 

Here we analyzed what happens when a cylinder or a 
disk rocks to an almost flat collision on its bottom sur- 
face. We found that the smallest deviation from a perfect 
face-down collision of a container's bottom results in a 
rapid rolling and/or sliding motion in which the contact 
point moves through a finite angle in infinitesimal time. 
Calculations of this finite angle explain certain apparent 
symmetry breaking in experiments involving rocking or 
toppled containers. 

In this system, the consequences of such a degenerate 
'collision' are a discontinuous dependence on initial con- 
ditions (rolling left or rolling right depend on the smallest 
deviations in the initial conditions). Such discontinuous 
dependence on initial conditions or geometry is a generic 
feature of systems in the neighborhood of simultaneous 
collisions. Other examples include a pool break or the 
rolling polygon of [H . 
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APPENDIX A: DERIVING EQUATIONS OF 
MOTION 

The moment of inertia tensor of the cylinder about its 
center of mass is, in dyadic form, 

I G = A(e x2 ® e x2 ) + A(e y2 ® e y2 ) + C(e z2 <8 e z2 ). (Al) 

The angular velocity SI of the cylinder is given by 

SI = ip sin e x2 — <p e y2 + (0 + ip cos 4>) e z2 ■ ( A2) 

The angular momentum about the center of mass is 

H/g = Ig SI 

= Aip sin (p e x2 — A<j) e y2 + C (6 + ip cos (f>) e z2 

(A3) 

The rate of change of angular momentum H iq is given 
by the sum of two terms: 1) the rate of change relative 
to the rotating frame e X 2-e y 2-e Z 2, obtained by simply 
differentiating the components in Eq. lA3[ and 2) the rate 
of change Sl 2 x H/ G due to the rotation of the e X 2-e V 2-e Z 2 
frame with an angular velocity Sl 2 = ip sin <f> e X 2 — <t> e y2 + 
-0cos0e x2 - 

H/g = (2Aip<j) cos <f> 4- Aip sin <p — Ccp{9 + ip cos <j)) \ e x2 
+ (^—Ap + Aip 2 sin 4> cos 4> — Cip sin </>(# + ip cos 4>)j e y2 
+ C (-^sin^ + i/Jcos^ + fl) e z2 (A4) 



1. No-slip rolling 

The point S on the cylinder in contact with the ground 
has zero velocity: vs = 0, no slip. Using rsc = rpc = 
Re X 2 + He Z 2, the velocity vg of the center of mass is 
given by 

vg = v s + SI x rpc 

= Q-H4>e x2 + R4>e z2 

+ [R9 + Rip cos 4> - Hip sin e y2 . (A5) 

The acceleration of the center of mass &g is given by 
adding two terms: 1) acceleration relative to the rotating 
e X 2-e y 2~e Z 2 frame, obtained by differentiating the com- 
ponents in Eq. IA5[ and 2) an acceleration term St x vq 
due to the rotation of e X 2-e y 2-e Z 2 frame. We obtain 

aG = dv G /dt 

= (-H4> - R<p 2 - Rip cos <p9 - Rip 2 cos 2 <p 
+Hip 2 cos (p sin 4>)e X 2 
+ ( —2Rip<j> sin <p> — 2Hip(f> cos <p + Rip cos <p 
—Hip sin <p + R9)e y2 
+ (Rp + Rip9 sin <p + Rip 2 sin <f> cos <j> 
-Hip 2 sin 2 (f>- H<f> 2 )e z2 . (A6) 
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The contact force Fp is given by linear momentum bal- 
ance: 

Fp = msLQ + mge z . (A7) 

Then, the moment M iq of all the external forces about 
the center of mass is 

M/ G = r G p x Fp = -r PG x F P 

= Bie x2 + B 2 e y2 + B 3 e z2 , (A8) 

in which, upon simplification, we have 

B\ = mH(—2Rip<psin(p — 2Hip<pcos<p 

+ Rip cos <p - Hip sin <p + R6) 

B 2 = m(- qH sin (p + H 2 4> + H Rip9 cos (p + 2HRip 2 cos 2 tj> 

- H 2 ip 2 cos (j) sin <f> + gi? cos + i? 2 + R 2 ip9 sin 

+ R 2 ip 2 sin<p cos (p - RHip 2 ) 

B 3 = -mR(-2R4>j>aiD.<f>-2Hij)j)C08<j> 

+ Rip cos (p - H ip sin <p + R9). 

Equating M/g an d H/g gives us the angular momentum 
balance equations of motion for no-slip rolling in Eqs. 11121 

2. Frictionless sliding 

No friction implies that the horizontal velocity of the 
center of mass G is a constant and can be set to zero 
by appropriate reference frame choice, without loss of 
generality. Further, noting that zq = R sin <f> + H cos cj>, 
we have, by differentiating twice: 

a G = z G e z 

= ( (R cos <p — H sin (p)(p~ (R sin (p + H cos (p)(p 2 ^ e z 

= sine/) (^(Rcosc/) — Hsm.(f))(f> — (R sin <j) + H cos (f))(j) 2 ^J e x2 

+ cos0 ^(i?cos0 — H sin <j))'(j) — {Rsincp + H cos c/))(j) 2 ^j e z2 

As before, we compute the contact force as Fp = mac + 
mge z and compute the net moment Mg about the center 
of mass as 

M/g = r G p x Fp 

= m(R cos (j) — H sin </>)(— R<f> 2 sin <j> 

+ Rep cos <j) — H(j) 2 costp — Hcpsinip + g)e y2 



Again, equating M/g with the H/g before gives us the 
required equations of motion Eqs. 1 1 13 1 



APPENDIX B: SCALING OF b 2 
1. No-slip rolling 

We now establish that b 2 = 0(<p 2 ), used in obtaining 
Eq. [151 The total mechanical energy E of the cylinder, 
a constant, is given by: 

2E = m(v G • v G ) + rj ■ Ig^ + P. E. (Bl) 
= m<f> 2 (H 2 + R 2 ) + C(9 + cos (pf 
+m(R8 + Rip cos <p — Hip sin (p) 2 
+A<P 2 + AiP 2 sin 2 <p + P.E. (B2) 

in which P.E. is the potential energy. We consider this 
total energy at the time of lowest tip angle <p, t = t c , 
we have <j> = 0, <p = <p c , ip — b 2 /(p 2 , and the potential 
energy P.E. = mg(Rsin <p+H cos (p) = Constant +0(<p c ). 
Making use of the usual small cp approximations (as used 
in the main text) and Eq. [7l namely ip + 9 ~ ip<f>, we find 
that 

S ~^=(|)%c- (B3) 

Because E is a constant, we have b 2 — 0((p 2 ) or b 2 = 
0(<p c ) as claimed earlier. 



2. Frictionless sliding 

The total mechanical energy E is now given by 

2E = m(p 2 (Rcos(p- Hsincp) 2 

+A1P 2 sin 2 <p + A<p 2 + C{6 + ip cos 0) 2 
+mg(R sin (p + H cos <p). (B4) 

Again, using small (p c approximations and considering 
the energy equation at t = t c , we obtain b 2 = 0((p c ), 
using arguments identical to the no-slip rolling case. 
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